home *** CD-ROM | disk | FTP | other *** search
/ Language/OS - Multiplatform Resource Library / LANGUAGE OS.iso / dsp / dr.bub / 96000.lha / 96000 / appb / b130.asm < prev    next >
Assembly Source File  |  1992-04-28  |  2KB  |  27 lines
  1. ; This program was originally published in the Motorola DSP96002 Users Manual
  2. ; and is provided under a DISCLAIMER OF WARRANTY available from Motorola DSP
  3. ; Operation, 6501 William Cannon Drive West, Austin, Texas 78735-8598.  For
  4. ; more information, refer to the DSP96002 Users Manual, Appendix B, DSP
  5. ; Benchmarks.
  6. ;
  7. ; B.1.30    Newton-Raphson Approximation for SQRT(x)  
  8. ;The approximation of sqrt(x) may be performed by using the  Newton-Raphson iteration to first find 
  9. ;1.0/sqrt(x).  The sqrt(x)  then can be approximated by x*(1.0/sqrt(x)).  
  10. ;
  11. ;
  12. ;              Newton-Raphson Approximation                  Program    ICycles 
  13. ;                     of SQRT(x)                               Words 
  14.  
  15.     seedr   d5,d4                    ;y approx 1/sqrt(x)    1       1 
  16.     fmpy.s  d4,d4,d2       #.5,d7.s  ;y*y, get .5           2       2 
  17.     fmpy.s  d5,d2,d2       #3.0,d3.s ;x*y*y, get 3.0        2       2 
  18.     fmpy    d4,d7,d2  fsub.s d2,d3  d3.s,d6.s ;y/2, 3-x*y*y 1       1 
  19.     fmpy.s  d2,d3,d4       d6.s,d3.s ;y/2*(3-x*y*y)         1       1 
  20.     fmpy.s  d4,d4,d2                 ;y*y                   1       1 
  21.     fmpy.s  d5,d2,d2                 ;x*y*y                 1       1 
  22.     fmpy    d4,d7,d2  fsub.s d2,d3  d3.s,d6.s ;y/2, 3-x*y*y 1       1 
  23.     fmpy.s  d2,d3,d4                 ;y/2*(3-x*y*y)         1       1 
  24.     fmpy.s  d5,d4,d4                 ;x*(1/sqrt(x))         1       1 
  25.  ;                                                          ---     --- 
  26.  ;                                                  Totals:  12      12 
  27.